Stochastic resonance in a simple model of magnetic reversals 
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We discuss the effect of stochastic resonance in a simple model of magnetic reversals. The model 
exhibits statistically stationary solutions and bimodal distribution of the large scale magnetic field. 
We observe a non trivial amplification of stochastic resonance induced by turbulent fluctuations, 
i.e. the amplitude of the external periodic perturbation needed for stochastic resonance to occur is 
much smaller than the one estimated by the equilibrium probability distribution of the unperturbed 
system. We argue that similar amplifications can be observed in many physical systems where 
turbulent fluctuations are needed to maintain large scale equilibria. 
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THE PROBLEM 

In this paper we discuss the effect of Stochastic Reso- 
nance (SR) for a simple model of magnetic reversals. As 
we will discuss later on, we show that the effect of SR 
can be amplified for systems where statistically station- 
ary states are observed, i.e. where metastable equilib- 
ria are due to non linear equilibration between turbulent 
fluctuations and non linear large scale effect. We argue 
that this effect is relevant in many physical systems and 
could be eventually observed experimentally. 

Before defining more precisely the problem we want to 
focus, it is worthwhile to review shortly the basic idea 
behind SR. SR was introduced almost 30 years ago in 
and within the framework of long term climate 
theory. One can be understood the mechanism of SR in 
the simple case of the stochastic differential equation: 



d(j) = (met) - g<j?)dt + y/adW{t) 



(1) 



where dW(t) is gaussian white noise 6— correlated in 
time. Because of the noise, <j) shows a bimodal probabil- 
ity distribution peaked around ±<fi m where (p^ = m/g. 
The average transition time t between the two peaks is 
proportional to 



exp ■ 



(2) 



It is well known that the transition time is a random vari- 
able exponentially distributed for small a. If we add on 
the r.h.s of ([1]) a periodic perturbation Asin(ujt), some- 
thing interesting can happen. For tt/lo ~ r the behavior 
of <j> becomes nearly periodic, i.e. 4> "jumps" between the 
two states ±<j) m periodically with period 27r/w. This be- 
havior can be understood in a number of different ways 
and we refer the reader to the original papers [l| and 
Q, see also 0] for a review on SR. For u> small and the 



deterministic time scale 1/m much shorter than r , the 
condition for SR to occur [lj can be written as 



A<j> Tl 



(3) 



In the simplified model ([T]), the equilibrium probability 
distribution P(<j>) is peaked around ±</> m which are stable 
stationary solutions of ([T]) for a = 0. In many physical 
systems, however, the probability distribution of the rel- 
evant order parameter (let us still call it ) is bimodal al- 
though there exists no stable stationary states: the peaks 
in the probability distribution arise because non trivial 
equilibration of internal dynamics which on the average 
can be approximated by an effective equation similar to 
([1]). This implies that the peaks in the bimodal proba- 
bility distribution should correspond to statistically sta- 
tionary solutions. A particular interesting case is the one 
where statistically stationary solutions are due to the bal- 
ance between non linear terms and internal fluctuations. 
In these cases, the effect of an external periodic perturba- 
tion can change the magnitude of the fluctuations which, 
in turn, change the parameters of the effective equation 
(PQ) i.e. the value of <fi m . Because of © , we can observe an 
amplification of SR. Such a mechanism is indeed observed 
in [4| for the case of two dimensiona Landau- Ginzburg 
equation. It is the purpose of this paper to describe the 
amplification of traditional SR in the case of a simplified 
model of magnetic reversals. Our major point is that, 
due to turbulent fluctuations in the presence of statisti- 
cally stationary solutions, SR is strongly amplified. Our 
example is just one of the many possible cases where 
the same amplification of SR may be observed and the 
present study outline the role of turbulent fluctuations in 
amplification of SR. We argue that other cases, relevant 
to turbulent flows and large scale dynamics of geophysical 
flows, may show similar amplification. 
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SIMPLE MODEL OF MAGNETIC REVERSALS 

The question of transitions between statistically so- 
lutions is central in the behavior of many out-of- 
equilibrium systems in physics and geophysics 0-11. As 
one particular example addressed here, we note that nat- 
ural dynamos are intrinsically dynamical. Complex mag- 
netic field evolutions have been reported for many sys- 
tems, including the Sun and the Earth Q. Formally, 
the coupled set of momentum and induction equations 
are invariant under the transform: (u, B) — ► (u, — B) 
so that states with opposite polarities can be generated 
from the same velocity field (u and B are respectively 
the velocity and magnetic fields). Such reversals have 
been observed recently in laboratory experiments using 
liquid metals, in arrangements where the dynamo cycle 
is either favored artificially fioj or stems entirely from 
the fluid motions [ll|, EH- I n these laboratory experi- 
ments, as also presumably in the Earth core, the ratio of 
the magnetic diffusivity to the viscosity of the fluid (mag- 
netic Prandtl number Pm) is quite small. As a result, the 
kinetic Reynolds number Ry of the flow is very high be- 
cause its magnetic Reynolds number Rm = RvPm needs 
to be large enough so that the stretching of magnetic 
fields lines balances the Joule dissipation. Hence, the dy- 
namo process develops over a turbulent background and 
in this context, it is often considered as a problem of 
'bifurcation in the presence of noise'. 

Building upon the above observations, we shortly re- 
view here a recent model proposed in fl3j ] which which in- 
corporates hydromagnetic turbulent fluctuations (as op- 
posed to 'noise') in a dynamo instability. We consider 
an 'energy cascade' model i.e. a shell model aimed at 
reproducing few of the relevant characteristic features of 
the statistical properties of the Navier-Stokes equations 
[li] ]. In a shell models, the basic variables describing the 
'velocity field' at scale r n = 2~ n r = fc" 1 , is a complex 
number u n satisfying a suitable set of non linear equa- 
tions (here ro = 2). There are many version of shell 
models which have been introduced in literature. Here 
we choose the one referred to as Sabra shell model. MHD 
shell model - introduced in [HI - allow a description 
of turbulence at low magnetic Prandtl number since the 
steps of both cascades can be freely adjusted [l||[I3|- We 
consider a formulation extended from the Sabra hydro- 
dynamic shell model: 

^ = l -(^ n (u,u)~^ n (B,B))- V k 2 n u n + f n , (4) 
at 3 

^ = |(*„(u, B) - $„(B, «)) - v m k 2 n B n , (5) 
where n = 1,2, ... and 

$ n (u,w) = fc n +i[(l + 6)u n+2 w* +1 + (2 - 5)u* n+l w n+2 ] 
+k n [(l - 28)u* n _ l w n+ i - (1 + <5)u n+ iu)*_ 1 ] 
+fc„_i[(2 - 6)u n ^iw n -2 + (1 - 26)u„-2W n -i] , (6) 



for which following [ll we chose 6 = —0.4. For this value 
of (5, the Sabra model is known to show statistical proper- 
ties (i.e. anomalous scaling) close to the ones observed in 
homogenous and isotropic turbulence. The model, with- 
out forcing and dissipation, conserve the kinetic energy 
Ey = S n |u„| 2 , the magnetic energy Eb = S„|B n | 2 and 
the cross-helicity i?e(S„M n B*). In the same limit, the 
model has a U(l) symmetry corresponding to a phase 
change cxp(i#) in both complex variables u n and B n . 
The quantity §> n (v,w) is the shell model version of the 
transport term vVw. The forcing term f n is given by 
fn = Si n fo/u*, i.e. we force with a constant power in- 
jection in the large scale. We want to introduce in eq. 
(J5]) an extra (large scale) term aimed at producing two 
statistically stationary equilibrium solutions for the mag- 
netic field. For this purpose, we add to the r.h.s. of (|SJ) 
an extra term M 2 (B 2 ), namely for n = 2 eq.© becomes: 

dB 

= F 2 (u,B) - M 2 (B 2 ) - v m k\B 2 (7) 

where F 2 (u,B) is a short hand notation for 
i/3($ 2 (u,B) - $ 2 (-B,u)). The term M 2 (B 2 ) is 
chosen with two requirements: 1) it must break the 
U(l) symmetry; 2) it must introduce a large scale 
dissipation needed to equilibrate the large scale mag- 
netic field production. There are many possible ways 
to satisfy these two requirements. Here we simply 
choose M 2 (B 2 ) = a m B\. From a physical point of view, 
symmetry breaking also occurs in real dynamos since the 
magnetic field is directed in one preferential direction 
which changes sign during a reversal. Also, large scale 
dissipation must be responsible of the equilibration 
mechanism of the large scale field. The choice of a 
non linear equilibration is made here to highlight the 
the existence of a non linear center manifold for the 
large scale dynamics. In other words, eq.(4) with 
M 2 (B 2 ) = a m B 2 is supposed to describe the 'normal 
form' dynamics of the large scale magnetic field. Note, 
that our assumption on M 2 does not necessarily imply 
a time scale separation between the characteristic time 
scale of B 2 and the magnetic turbulent field. Finally, 
since the system has an inverse cascade of helicity, we 
set B\ = as boundary condition at large scale in order 
to prevent non stationary behavior. 

The free parameters of the model are the power input 
/o, the magnetic viscosity v m and the saturation parame- 
ters a m . Our numerical simulations have been done with 
n = 1, 2.., 25. Actually, the parameter /o could be elim- 
inated by a suitable rescaling of the velocity field. We 
shall keep it fixed to /o = 1 — i. In this system, a possi- 
ble estimate of Reynolds numbers is Ry = J (Ey) jk 2 v = 
^(£V)r /4zy and R M = y/(E v )r /4u m . 

For very large v m , the magnetic field does not grow. 
Then, for v m lower than some critical value, (B 2 ) as well 
as Eb increases for decreasing v m . Eventually, ( \B 2 1) sat- 
urates at a given value while Eb still increases, showing 
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FIG. 1: Time behavior of i?2r for three different values of v m 
(displayed on the left side) and constant v. The blue segment 
in the upper panel shows lOIMd, where td is the dissipative 
time scale computed as td = l./(k^v m ). One time unit in 
the figure corresponds to the large scale eddy turnover time 
Numerical simulations have run for much longer 
than the time intervals shown here - in the complete time 
series there is no asymmetry in between the ±i?2 states. 



the upper panel of figure [T] is much longer than the char- 
acteristic time scale of B2 near one of the two equilibrium 
states. The system spontaneously develops a significant 
time scale separation, for which given polarity is main- 
tained for times much longer than the magnetic diffusion 
time, fn figure [5] we show the average reversal time as a 
function of v m . 
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that for v m small enough a fully developed spectrum of 
B n is achieved. This type of behavior is i n agree ment 
with previous studies of Taylor- Green flows 
flows in a sphere [2l[ or MHD shell models 
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The onset of dynamo implies that there exists a net flux 
of energy from the velocity field to the magnetic field. At 
the largest scale, the magnetic field B2 is forced by the 
velocity field due to the terms Fi{u,B). The quantity 
S = 1l\F2{u, B)B%\ in (fT0|) is the energy pumping due 
to the velocity field which is independent on B2 and a m . 
Thus, from eg. (|13|) we can obtain: 

ld\B 2 



2 dt 



= S- a m \B 2 \ 2 (B 2 2r - B 2 2l ) - is m k 2 2 \B 2 \ 2 (8) 



where i?2r and B2i are the real and imaginary part of 
B 2 . For large u m , the amplitude of B2 is small and the 
symmetry breaking term proportional to a m is negligible. 
Under this condition, and with the boundary condition 
constrains, we expect from © or (fTU|) that the behavior 
of B 2 is periodic, as it has been observed in the numerical 
simulations. On the other hand for relatively small v mi 
the non linear equilibration breaks the U(l) symmetry 
and B 2 i becomes rather small and statistically station- 
ary solutions can be observed with B\ r = y/ S/a m . In 
[131 ] . it is discussed a systematic study of the magnetic 
reversal as a function of v m . In figure [U we show three 
different time series of the i?2r = Re(B 2 ) as a function 
of time for three different, relatively large, values of the 
magnetic diffusivity. The figure highlights the two major 
informations, namely the obervation of reversals between 
the two possible large scale equilibria and the dramatic 
increase of the time delay between reversals for increasing 
v m values. Note that this long time scale, as observed in 



FIG. 2: Average persistence time r as a function of the mag- 
netic viscosity i>m for a m = 0.1 and v = 1CF 7 . The green line 
corresponds to the fit given by equation (|12[1 . In the insert 
we plot l/ln(r) and its error bars versus v m to highlight the 
linear behavior predicted by (|12[) . Note that the error bars 
are smaller than the symbol size except for the very last point 

In order to develop a theoretical framework aimed at 
understanding the result shown in figure [H we assume, 
in the region where (|i?2| 2 ) is independent on v m , that 
B 2 i ~ and that the term F 2 (u,B) can be divided into 
an average forcing term proportional to i?2r and a fluc- 
tuating part: 



F 2 (u,B) = PB 2 



(9) 



where (3 depends on /o and 4>' is supposed to be uncor- 
rected with the dynamics of B2 , i.e. ([</>' -BJ]) = 0. Note 
that in the context of the mean- field approach to MHD, 
the first term /3i?2 would correspond to an 'alpha-effect'. 
Using ([5]), we can rewrite the equations for B2 as follows: 



dB 2 
~~dt 



(10) 



where we neglect the dissipative term since /3 3> v m k\ 
in the region of interest. eq. fTU)) must be considered an 
effective equation describing the dynamics of the mag- 
netic field B 2 and its reversals, and the fluctuations <// 
incorporates the turbulent fluctuations from the velocity 
and magnetic field turbulent cascades. It is the effect 
of (j)' which makes the system 'jump' between the two 
statistically stationary states. Using (jHJ we can obtain 
P = y/Sa m while the two statistical stationary states can 
be estimated as ±Bq, Bq = /3/a m . 
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FIG. 3: Solutions of eq. (|13l as compared to the behavior of 
B 2r obtained from Q-JSll for v m = 0.00026. The choice of 
a in (|13[) is chosen to reproduce the mean transition time r 
close to the solution of j4])-([5]) 

Interpreting cq. (fTU|) as an effective stochastic differen- 
tial equation, we can predict r to be 







= exp (|j 


\a m <? J 





where a is the variance of the noise 4>' acting on the 
system. In [l3| it is suggested that a = A(v^ — v m )/uL 
which leads to 

where C is a constant independent of v m . This func- 
tional form is displayed in figure [2j it agrees remarkably 
with the observed numerical values of r for a rather large 
range. In the insert of figure [T] we show l/log(r) as a 
function of v m to highlight the linear behavior predicted 
by eq.(fl"2|). The physical statement represented by (fT2|) 
is that the average reversal time should show a critical 
slowing down for relatively large v m . In other words, we 
expect that fluctuations around the statistical equilibria 
increase as Rm increases. The increase of fluctuations 
may not be monotonic for very large Rm , which explains 
why we are not able to fit the entire range of v m shown 
in figure (Q]). 

THE EFFECT OF SMALL PERIODIC FORCING 

As shown in the previous section, eq. (fTOjl can be con- 
sidered an effective stochastic differential equation. More 
precisely, we can approximate the reversals of magnetic 
field Bi T by using eq. ([T|) with (f> = B-2 r and suitable val- 
ues of m and g such that <j) m corresponds to the magni- 
tude of the observed statistically stationary states in the 
probability distribution of i?2r- Finally we must choose 
(7 in ([I} such that the average transition time is equal 



FIG. 4: Numerical solutions of (|14[) for different values of 
the external periodic forcing Asin(2nt / u) . For A — 0.4 the 
dynamics of <f> starts to behave periodically 

to the one observed for the magnetic reversals (see fig- 
ure Q). For our purpose, we choose v m = 0.00026 so 
that the average transition time is order 4000. Then the 
effective equation ([TJ is 

d<f> = (j){(j) 2 m - <f> 2 )dt + ^dW{t) (13) 

where 4> m = 1.51 and a = 0.34 are the numerical values 
chosen in (| 13[) for the dynamics of <j> to be close to the 
observed numerical behavior of i?2r- In figure we 
show the behavior of the numerical solution of (|T3l) as 
compared to the time behavior of Bir computed for v m = 
0.00026. Hereafter we refer to eq. (fl~3|) as "double well" 
model while we use the term "mhd" for eq.s (jj)-©. 

Following our discussion in the introduction, we can 
use ([3|) to estimate the amplitude of an external periodic 
forcing Asin(u)t) applied on the r.h.s of (|13|) with 2it/uj = 
4000. It turns out that © gives A ~ 0.5 for SR to occur. 
This agrees very well with the numerical results shown 
in figure (0| where we plot the numerical solution of 

d4> = (j)((t) 2 n - (j) 2 )dt + v^dW(t) + Asin(2n/ujt) (14) 

with A = 0.1 (upper panel) A = 0.2 (middle panel) and 
A = 0.4 (bottom panel). A more quantitative descrip- 
tion can be obtained by looking at the average Fourier 
amplitude Pa = (|<Mw)|) and plotting Pa as a function 
of A. This is done in figure © (red circles). 

We now turn our attention to the system @ -([5]) which 
describes the magnetic reversals within the simplified 
model discussed in the previous section. We add to the 
equation of motion of Bi r an external periodic forcing 
Asin{2irt / lo) with the same period 2-Kt/uj = 4000 chosen 
for the numerical simulations of figure ((3|). In figure ([5]) 
we show the behavior of Bi r as a function of time for 
A = 0.08 and compared it against the solution of (TT41 for 
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FIG. 5: Upper panel: numerical behavior of Bi r obtained 
from (|4j|-([5jl when an external forcing Asin(2nt / ui) is added 
to r.h.s. of ((5} for n = 2. The blue curve shows the periodic 
forcing. Lower panel: solution of (|14|l with A — 0.1. 
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FIG. 6: The average Fourier amplitude Pa = {\4>(co)) as a 
function of A for eq. (|14|) (red circles). The blue triangles 
refer to the same quantity computed for the eq.s - d5j) : the 
response function is much larger than for the case of stochastic 
differential equation (|14[) . 

the same A. We can clearly observe SR in the mhd equa- 
tions (|4])-[5]) while for eci.(fT4f the effect of A is too small. 
In figure ([5]) we show the the average Fourier amplitude 
Pa = (|i?2r(w)|) and plotting Pa as a function of A com- 
paring the result with the same quantity computed for 
the stochastic differential equation (jT4")) . The conclusion 
is that SR is amplified in the mhd equations (HI)-© by 
almost a factor 5! 

We remark that the effect shown in figure ([5]) is higly 
non trivial, i.e. it is not trivial to figure out an effec- 
tive stochastic differential equation similar to (p~3|) which 
shows the same sensitivity to the external perturbation 
A. Using a different language, borrowed by statistical 
field theory, we can say the the effective eq. (|13[) is de- 
scribed in terms of renormalizcd parameters which de- 
pends on the turbulent fluctuations. Since the statisti- 



cally stationary states are due (on the average) to non 
linear equilibration between the fluctuating forcing term 
F2(u,B) and the large scale term a m £?2, the effect of 
external perturbation changes (probably in a non lin- 
ear way) the amount of fluctuations which fix the values 
of the renormalized parameter in (|13|) . Another impor- 
tant point to remark is that the behavior observed in 
our simplified mhd model can be blindly parameterized 
as "stochastic nose", i.e. although turbulence is quali- 
tatively acting as a noise in the dynamics, the turbulent 
fluctuations are correlated to the statistical equilibria in a 
way which is hardly to parameterize as an external noise. 

CONCLUSION 

In this paper we have shown a relatively simple ex- 
ample of amplification of SR in a system characterized 
by statistically stationary states. It is interesting to re- 
mark that, in our system, the large scale magnetic field 
is fluctuating in time around some state which is fully 
maintained by the turbulent energy flux from the velocity 
field (i.e. the dynamo instability). The effect on exter- 
nal perturbation changes this equilibria and amplifies the 
susceptibility of the system. The amplification observed 
shows that turbulent fluctuations cannot be parameter- 
ized as "noise" independent of statistical equlibria. 

We argue that there may exist many different phys- 
ical systems where similar effects can be observed. In 
particular, it will be interesting to explore whether such 
a strong sensitivity to external perturbation is relevant 
in geophysical flows where many theories of large scale 
multiple equilibria have been proposed in the past. Also, 
concerning the climate theory, our simple but non triv- 
ial example, shows how non linear fluctuations could be 
coupled to external forcing in a rather non intuitive way. 
This may open the possibility to reconsider SR in cli- 
mate theory in the framework of more complex models 
where internal turbulent fluctuations of climate dynamics 
are explicitly taken into account and simulated. Finally, 
within the application of our present result in the case 
of dynamo instability, we argue whether a similar effect 
could be observed experimentally, i.e. whether by apply- 
ing a small external forcing a large SR is observed for 
relatively low value of the external amplitude. 
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